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' Abstract 

O ' 

Time-dependent Hartree-Fock simulations of the evolution of ex- 
cited gold fragments have been performed. The observed dynamics 
appears more complex than the collective expansion picture. The min- 
imum density is often not reached during the first density oscillation 
because of the beating of several collective compression modes. 



1 Introduction 

Studying the dynamics of a highly excited nucleus, such as those created 
in Heavy-ion collisions, is a priori a hard task since we are dealing with a 
many-body quantal problem involving a large number of degrees of freedom. 
A widely used approximation is to reduce the problem to the dynamics of 
few collective variables, the remaining degrees of freedom being treated in a 
statistical manner. In particular, excited nuclear systems are often assumed 
to simply be hot nuclei expanding or contracting isentropically in a unique 
collective compression mode, the so-called breathing mode ( see for example 
the expanding fireball scenario or the fragmentation models [ffl 0, [3]]). 

This picture has been recently applied to explain the GSI data concerning 
the correlation between the apparent temperature and the excitation energy 
of fragmentation events known as the Aladin "caloric-curve" (see ref.[[| and 
refs. therein) . In the approach of ref . M the fragmentation of various nuclear 
systems formed by the projectile remains after a collision, is assumed to occur 
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at the respective turning points of the induced giant breathing mode. The 
location of these turning points in temperature and density is proposed to 
be the explanation of the observed " caloric-curve" || . More generally, most 
of the multifragmentation data are understood assuming that the system 
expands to reach a low density value either in statistical models || or in 
dynamical simulations [[F| using stochastic approaches ||, The general belief 
is that the system breaks during the first expansion associated with the giant 
breathing mode of the hot nucleus. More generally, this assumption is often 
made in simple semi-classical simulations or in multifragmentation models 

BBS- 

However, recent results about photon production during Heavy-Ion col- 
lisions seems not to be compatible with this simple scenario of a fast multi- 
fragmentation during the first expansion of the system and the data have been 
understood assuming at least one recompression before break-up j|T0|, ITT . 



On the other hand, from the study of monopole vibrations, it is well- 



known that nuclei can exhibit a large variety of compression modes [12], [13 1 



Moreover, microscopic description (using TDHF or is linear version) of monopole 
vibrations predicts that these modes cannot be simply described by the 
breathing mode picture. Since these giant resonances are the small ampli- 
tude limit of the collective expansion assumed to occur in heavy ion collisions, 
one may worry if the picture of a unique large amplitude collective breathing 
mode leading the hot nuclei to low density regions is not too simple. There- 
fore, in this article, we propose to study the large amplitude motion of an 
excited nucleus within the framework of the Time-Dependent Hartree-Fock 
(TDHF) theory in order to critically discuss the validity of the single breath- 
ing mode picture. It should be noticed that, heavy ion collisions are rather 
well described by simulations using the semi-classical version of TDHF: the 
Vlasov equation, possibly extended to take into account the damping com- 
ing from nucleon-nucleon collisions (i.e. the so-called BUU approach [|i4f ). In 
these simulations the expansion phase appears to be dominated by the mean- 
field dynamics confirming that TDHF can be a good tool to investigate this 
large amplitude motion. Furthermore, as far as we know, TDHF is the only 
well-tested quantum model tractable for Heavy-Ion collisions. Of course, 
a quantum model that include collision effects (extended TDHF) would be 
preferable. However, we believe that the introduction of this effect will not 
change considerably the discussion since, during the first expansion, the gain 
and lost term of two-body collisions are almost equal and so cancel. 

In this article, we will analyzed the TDHF dynamics and we will point 
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out that some aspects, important for the dynamics, cannot be accounted for 
in the simple breathing mode picture. This may change our understanding 
of the expansion dynamics of nuclei created during heavy ion collisions. In 
particular, we observe that, in general, the minimum central density is not 
reached during the first density oscillation. This may allow the system to un- 
dergo several monopolar oscillations before breaking in pieces as the photon 
data are suggesting. 



2 Mean-Field dynamics 

In mean-field evolution all the information is contained in the single-particle 
density operator p. Let us introduce the single-particle basis which diago- 
nalizes this operator: p = YJ 1^/3 > np < <&p\, where (3 contains all quantum 



numbers. Following the work of ref. [15], we consider the mean-field evolution 



of a spherical symmetric^], spin-isospin saturated nucleus. The mean-field is 
parameterized as 



U(r) = \t Q p{r) + ^±^p(rf +1 + V J ^z^p(r')d r> + V c (1) 

a 

The two first terms are the same as in ref. JTB ] and a general discussion could 
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be found in [16|. The Yukawa potential is introduced in order to mimic the 



gradient part of the potential as in ref. |17], 18| . Finally, Vc is the direct 



Coulomb potential which is ,in our calculation, approximately treated by 
given a charge to each nucleon. Parameters are chosen in order to have the 
same property as the SkM* [|19| force for the nuclear matter. We have taken 
t = -2191.73 MeV fm 3 , t 3 = 18818.8 MeV fm 7 / 2 , a = §, a = 0.45979 fm, 
Vq = —461.07 MeV fm 3 . For the nuclear matter, this leads to E/A = —15.8 
MeV at a saturation density = 0.16 fm -3 and to an incompressibility 
modulus of Koo = 200 MeV which corresponds to a soft Equation Of State. 
It should be noticed that the presented results do not qualitatively depend 
upon the compressibility of the nuclear force. 

1 Either experimental evidence, for example in central events selected with INDRA, 
or theoretical simulation of the BUU type confirms the fact that a spherical equilibrated 
source can be a reasonable ansatz. 
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Since in TDHF, the initial symmetry of the system is conserved, the 
dynamical evolution reduces to a set of equations for the radial part (p n i(r, t) 
of the single-particle wave functions <&p(r, t) 

ih-vrair) = — - — + y 2 ; + U(p(r, t)) Vnl {r) (2) 
at [ 2m or 1 2mr z J 

where n and I are the principal and orbital numbers and the density takes 
the particular form 

p(r,t)=4E(2l + lK i ^ 1 (3) 

n,l 

Numerically, we consider 40 orbitals, (p n i. The Schrodinger equation is 
solved in coordinate representation. The lattice of size 300 fm is discretized 
in steps of size Ar = 0.2 fm. An imaginary time method is used to generate 
the ground state of the static problem |L5L |l6f ( a small external field Ar 2 with 



A = 0.25 MeV fm -2 is added during this initialization in order to define the 
particle states). Occupation numbers are calculated afterwards according to 
the Fermi-Dirac distribution at the considered temperature. We perform the 
TDHF evolution up to 1500 fm/c with a time step of 0.75 fm/c. 

In order to study nuclei comparable to those expected to be produced in 
reactions, we take the same initial conditions as in ref. for the fragments of 
a Au-projectile. Let us first focus on one typical initial condition : an excited 
fragment of mass A = 191. The excitation energy is taken from the abrasion- 
ablation model[]20[]. Here we consider that the excitation energy is related to 
the number (AA) of particles lost during the abrasion-ablation stage of the 
collision by the relation E* = 13.3AA MeV. Furthermore, following ref. 0] 
we assume that the projectile fragments are initially slightly dilated. We 
implement this dilution by rescaling the wave-function according to r — > 
r 0.8 -1 / 3 . In this condition the temperature is adjusted in order to obtain 
the correct excitation energy. For the considered fragment which is closed 
from the projectile, the excitation energy is small and partly stored in the 
initial dilution so that the temperature is only T = 0.55MeV. As we will 
discuss in the following, the conclusion we will draw are rather general and 
the complex dynamics is observed for a wide range of initial value of the 
mass, the density and the temperature. 
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3 Results 



3.1 Density profiles 

In we display the density (solid line) at various times as predicted by 
the mean-field approximation^. We observe that the density profile display a 
lot of distortions. To emphasize this point we have fitted the density profile 
at each time by a Fermi shape 

P F (r,t)= P{t L R{t)) (4) 

Typical results of the fit are displayed in fig.[l](dashed curves). At initial time, 
the density profile is well fitted by a Fermi shape. In a scaling picture of the 
evolution, one expects, at all time, that this fit should be also good. The 
strong difference between the fitted Fermi shape and the real density profile 
at time t = 108 fm/c demonstrates how far we are from a simple scaling 
picture of the expansion. 

From the three parameters p(t), R(t) and the diffuseness coefficient a(t), 
we can define three scaling parameters as 

R(t) a(t) . ,.,.„, , 

R(t ) a(t ) \ p(t) 

where to is the initial time. If the collective dynamics is dominated by a 
unique breathing mode (a global scaling of the density profile) the three a 
should be equal and should oscillate with a unique frequency close to the one 
of the monopolar vibration of the considered nucleus.^ 

In fig.|] (Top), we show the evolution of the a's as a function of time. 
From this figure, it is clear that only a R and a p display regular oscillation 
with a given frequency. Performing a Fourier transform of this quantities, 
we found a marked peak at Hu = 12.5 MeV (which is around the expected 
breathing mode frequency of the considered nucleus at zero temperature). 

2 It should be noticed that we have checked that putting the ingredients considered in 
reference [13] we recover their results. The difference is that we are making a more detailed 
analysis for a heavier nucleus. 

3 A possible shift in frequency can be due to the anharmonicity of the large amplitude 
motion and to the effect of the temperature. 
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On contrary, a a contains clearly a complex superposition of many modes 
and is never comparable to the two other scaling factors. The Fourier trans- 
form of a a presents not only the frequency hu = 12.5 MeV but also two 
other broad peaks around 17 MeV and 26 MeV, and a small peak around 
40 MeV. The Fermi shape fit being unrealistic for the density profile at time 
t = 108/m/c, it is difficult to really conclude on the link between the pres- 
ence of different collective modes and the appearance of hollow shapes. In 
the following, we will characterize more precisely this link. 

From fig. p] it is clear that not only the surface of the nucleus is following a 
complex dynamics but also the interior which exhibits a tendency to produce 
hollow structures. To quantify this feature, we have computed two different 
averaged densities as a function of time: the solid line in figure ^] (bottom) 
represents the averaged density in a sphere of 2fm radius (p 2/m ) whereas the 
dashed line stands for the density averaged in a sphere of 5fm (p ).The 
sphere of radius 5fm is big enough to contain most of the volume of the 
nucleus and small enough (see figure |]) to not be polluted by the surface 
shape. As expected, p presents a regular oscillation at the breathing 
mode frequency. On contrary, it is clear that p 2/m contains more than the 
breathing mode, in particular, this quantity is more sensitive to the presence 
of hole at the center of the nucleus as seen in fig.0. Therefore, it appears 
that in addition to the expected large amplitude volume mode other modes 
affecting both the interior and the surface of the nucleus are present in a 
complete TDHF simulation. 



3.2 Collective vibrations 

In order to quantitatively analyze the observed complex oscillation pattern of 
the density, we have performed the Fourier transform of p(r, t) leading to the 
spectral density p(r, u). The Fourier transform is performed up to 1500/m/c 
which is sufficient to have a resolution of few hundred keV. A plot of |p(r, u>)\ 
for p/p = 0.8 is displayed in fig. |3|(Top). We see in particular that non-local 
waves develop during the evolution. Of course we recover the typical breath- 
ing mode at 12.5 MeV but also contributions around 17MeV and 26MeV. 
The radial dependence of p(r, uj) for a given u can be considered as the tran- 
sition density of the corresponding mode. Therefore, the appearance of hole 
at the center of the nucleus seems to be connected with the surface waves. 
We have tested that removing the coulomb field does not modify the pre- 



sented conclusions (as suggested in [PqQ. Therefore, the creation of a hollow 
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structure is not a consequence of the coulomb interaction. 

It should be noticed that starting from a value of 0.8, p 2fm /p reaches 
much lower values during its evolution. Moreover, the minimum minimorum 
(namely p 2fm /p = 0.4) is not reached during the first oscillation but at 
the third ones. The observed oscillation pattern of p 2fm is characteristic of 
the beating of several modes as clearly seen from the Fourier transform. It 
should be noticed that reducing the amplitude of the oscillations we have not 
observed a major modification of the Fourier spectrum. An example of such 
a case is shown in figure 3 (Bottom) for initial compression p/p = 1. As 
far as the main peaks are concerned, the frequency structure of the Fourier 
transform does not change by more than a half MeV when the amplitude of 
the mode increased. However, looking in details we can see that the relative 
importance of the different modes is not the same. This is what we expect 
since the amplitude of each mode is directly related to the initialization 
procedure. This demonstrates that the coupling between modes and the 
anharmonicities are small, and that the main effect is the beatings of different 
modes. However, looking to the details of the radial structure of the Fourier 
transforms we can notice that the transition densities are slightly modified by 
the modification of the initial amplitude. This is clearly evidenced in figure 
4, where the transition densities p(r, u) for the two main peaks u = 12 MeV 
and uj = 17 MeV is displayed for the two initializations shown in figure 3. In 
this figure we can see that the Giant monopole resonance around 12 MeV is 
more robust than the smaller peak at 17 MeV for which a contribution in the 
interior of the nucleus is seen at large amplitude. This volume contribution 
explains why a hollowed shape can appear as a function of time as the two 
modes are beating one again the other with different nodal structures. 

4 Conclusion 

In the presented simulation, the beatings of different monopole modes lead 
to very low densities in the interior of the nucleus and to the formation of a 
hollow structure after a long time (2 or 3 oscillation period). This tendency 
of excited nuclei to expand in hollow structures have already been observed 
in several dynamical calculations of the Vlasov type |22|, |23| however the new 
feature shown by the presented quantum calculations is that this structure 
gets more pronounced after several collective oscillations. This analysis shows 
that this is due to the presence of several collective degrees of freedom with 
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a transition density affecting the interior of the nucleus. 

This could change the discussion about the possible fragmentation of nu- 
clei and also about the involved time scales. Indeed, a possible interpretation 
of multifragmentation is that heavy excited system could reach dynamically 
a configuration in which the nucleus is unstable against a splitting in pieces. 
This might be in particular the case for low density regions or for hollow 
structures. In the framework of one collective variable picture, after pre- 
equilibrium, the excited nucleus reaches the lowest density region at the first 
turning-point of the collective expansion. This has given a widely admitted 
time-scale of this phenomenon and is now often used to discuss break-up 
of a systemQ. If now we focus on P2f m , we see that a nucleus could oscillate 
several times before reaching a configuration presenting density profiles with 
low density regions. If this configuration is unstable we see that the nucleus 
could wait a long time, of the order of 100 fm/c for the second turning-point 
or even 200 fm/c for the third, before initiating its fragmentation. It should 
be noticed that indications in direction of this scenario have been recently 
reported in hard photons study of Heavy- Ion collision [IT]. 

Note also that we have observed the same behavior for different nuclei 
with higher temperatures. Moreover, we have applied the same analysis with 
a hard EOS (SIII force), and the same observations have been made. There- 
fore, the reported conclusion are rather generic. From the analysis we have 
related this behavior to the existence of a volume contribution with several 
nodes in the transition density. This, together with the Landau spreading, 
may occur as generic features and explains that the presented dynamics was 
found in many simulations. 

In summary, we show that realistic TDHF calculations of expanding nu- 
clei cannot be reduced to the simple large amplitude breathing mode picture. 
Generically, the observed dynamics presents the beating of several modes so 
that it may happen that low density regions can be reached after several 
global oscillations. This may have strong qualitative and quantitative conse- 
quences on the dynamical expansion of nuclei: in particular, time-scale could 
be strongly increased and exotic shapes can be reached. 
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Figure 1: Density profile (solid curve) for different time. Dashed curve rep- 
resent fits with a Fermi shape parametrisation of the density. This fit is 
not able to reproduce the complicated density profile (see t = 108 fm/c) 
and demonstrate the difference between the TDHF evolution and the scaling 
picture. 
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Figure 2: Top: Evolution of the three scaling factor calculated with the three 
parameters of the Fermi shape density. Dashed line: a a , long-dashed line: 
qir and solid line: a p . Bot: Evolution of the central densitie P2f m (solid 
curve) and p5/ m (dashed curve) as a function of time. 
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Figure 3: Fourier transform \p(r, ou)\ (3D view). Two initial conditions are 
studied p/po = 0.8 (top) and p/p — 1 (bottom) are drawn. Fourier trans- 
forms are performed over 1500/m/c. 
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Figure 4: Transition densities p(r, a;) for u = 12 MeV (solid line) and uj = 17 
MeV (Dashed line) drawn for two values of the initial compression (see fig 
3). 
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